Exponential spreading and singular behavior of quantum dynamics 

near hyperbolic points 
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Quantum dynamics of a particle in the vicinity of a hyperbolic point is considered. Expectation 
values of dynamical variables are calculated, and the singular behavior is analyzed. Exponentially 
fast extension of quantum dynamics is obtained, and conditions for this realization are analyzed. 
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Hyperbolic (saddle) points are a source of instability 
in dynamical systems fl)]. Therefore, quantum dynam- 
ics of a particle in a saddle point potential is an impor- 
tant problem in quantum chaos [2tJ1- The hyperbolic 
point at the origin {x,p) = (0,0) can be described lo- 
cally by the Hamiltonian Hioc = xp. The Lyapunov ex- 
ponents detecting stable A_ and unstable A+ manifolds 
are A± = ±1. This leads to the exponential spreading in 
quantum dynamics and exponential growth of observable 
quantities [3, Q ■ The Hamiltonian Hioc has been studied 
in connection with the Riemann hypothesis 0,[1|j scatter- 
ing of the inverted harmonic oscillator [9| , and eigenstates 
near a hyperbolic point [10]. 

We consider a quantum system which is a second order 
polynomial of Hioc with the Hamiltonian 



H. 
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where w and fi are the linearity and nonlinearity pa- 
rameters correspondingly, while the coordinate and mo- 
mentum operators obey the standard commutation rule 
[x,p\ = ih with Planck's constant h. This system was 
considered in Q in connection with the Riemann hy- 
pothesis as well. Our primary interest in this Hamil- 
tonian is related to a problem of quantum dynamics con- 
sidered in the Heisenberg picture [llj where the expec- 
tation values of the operator s yi i) and x{t) were calcu- 
lated in the coherent states [12|, specially prepared at 
the initial moment t = Q. As shown in pjj, the dy- 
namics of these expectation values becomes singular at 
specific singularity times ti. For example, for the ob- 
servable value of x^ these singularities occur at times 
^' ~ ^ ^T^WTl^ ' ~ ■ • ■■ A remarkable prop- 



erty of these explosions is a pure quantum nature: in 
the classical counterpart this corresponds to the separa- 
trix motion without any singularities. As follows from 
the analysis of Ref. , this explosion behavior results 
from the interplay between the nonlinear term and the 
specific choice of zero boundary conditions on infinities. 
When /i = 0, the singularities are shifted to infinity and 
expectation values of operators are finite, and this result 
is independent of the boundary conditions. 

We develop a different consideration of the problem to 
understand the nature of this singularities that, as will be 
shown, are related to the choice of the initial conditions. 



First we consider the quantum dynamics of a particle 
of a unit mass in the saddle potential described by the 
Hamiltonian Q. For calculation of the expectation val- 
ues, following [ill, the initial wave function is chosen in 
the form of the coherent state ^'o(a;) = (xjo). In the 
x-representation this is the Gaussian packet (ill . [l2| : 

Vl'o(x) = {x\a) = (?,^)-l/4e-|-lV2fi-(.^-2y2.a+a^)/2;.^ 

(2) 

where a € C. It is worth mentioning that in this notation 
the dimension of x is y/h. For simplicity we calculate the 
expectation value of the operator x^{t). Thus we have 



{x\t)) 



%[x)U\t)x^U{t)'^o{x)dx, (3) 



where U{t) is the evolution operator. 

The axis of integration we split into three intervals 
(— oo,— cco], [— xo,a;o], [xq, oo), and the expectation value 
is expressed by the following three integrals: 

{x^t)) = hit) + ij{t) + i+{t) = r-i, ^i{t)x^-^s{t)dx 

+ /r^° ^){t)x^-^f{i)dx + Q ^}{t)x^^f{t)dx . (4) 

The dynamics in the finite interval [— xo, xq\ is considered 
in the framework of the truncated interaction. Near the 
hyperbolic point this dynamics is considered locally, such 
that H = Hs for x < xq and the particle is, for example, 
free with H = Hf outside the interaction region |a;| > xq, 
where Xq > determines arbitrary the interaction range. 
Note that we do not consider a scattering task and just 
truncate the integration. Here Hg is determined by Eq. 
([T]), while Hf = p^ 12 determines free motion. Therefore, 
the dynamics of an initial wave function is determined 
in these two different regions 



*/W = C^/W*o and *«(t) =t/,(t)*0: 
where the evolution operators Uf(t) = exp 



(5) 
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and Us{t) — exp —j-Hgt describe two independent pro- 
cesses, and corresponding shift of the wave functions is 
supposed. 

Integrals are calculated by substituting Eqs. (O in 
Eq. Q and taking into account the explicit form of the 



2 



Hamiltonians. In the free motion window the evolution 
of the square coordinate operator is 



u}{t)xUfit) 
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Therefore, the free motion integrals (t) do not have any 
particular features. Their values can be expressed in the 
form of the error function $(z) — (27r)~^/^ ^ e'' /^dyy 
. Then we obtain that (t) ^ , as expected. 
The saddle point integral possesses a more interest- 
ing behavior. Using calculations performed in [Eg. 
(4.16)], we arrived at the following expression: 



X exp 
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2ft / dx 

First, we admit the exponential quantum growth, ob- 
tained in Ref. [llj . This quantum behavior has clas- 
sical nature of the near separatrix motion, observed 
also for a kicked system At the singularity times 
~ 35H^ expression (O is simplified, and the 

integral is calculated exactly. It reads 
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where — a* +i{~iya. This expression is finite for finite 
Xq. When xq approaches infinity the integral diverges 
and ti are the singularity points. Note that xq is an 
arbitrary defined scale [14| . 

To generalize the consideration of the explosion sin- 
gularities, first we consider eigenvalue problem for the 
Hamiltonian Hs- Since the operator xp — ih/2 commutes 
with Hs, this problem is reduced to the dimensionless 
equation for the eigenfunctions Xei^)- 



with the following solution 



VNs 



: exp[ieln |x|] , 



(9) 



(10) 



which satisfies the boundary conditions Xe{x = ±c») = 
and Ns = Airh}/'^ TtI]. For the continuous spectrum the 
normalization condition is 



Xl'{x)xe{x)dx = (5(e - e') 



(11) 



(see e.g. |T5|). 

Now, expanding ^'o(x) over the complete set of x^i.^), 
we obtain 



*o(a:) 



deq{e)Xe{x) . 



(12) 



Note, that the explicit form of the expansion coefficients 
q{t) is not important, since integration over energy e will 
be performed with exactly the same form of Xc {x) . Sub- 
stituting Eq. (dH) in the integral Ig in Eq. ^ with 
Xq = oo, we obtain 



If{t) =21 X 
'o 

X X*e'{x)Xe{x)dx , 



q* {e')q{e)e~i{E - E')tdede' 
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where E = 2uje — Ahfit^ is the energy of Hg , and we use 
that Xe(~2;) = x<i{x)- The complex Gaussian exponents 
are presented in the form of the Fourier integrals: 



e^^^dT 



exp{=FjTVl6M} • (14) 



Substituting these expressions in Eq. (jl3p and taking 
into account the explicit form of Xeix), we obtain In |x| — 
2ujt — T = ln(|x|e~^"*e~'^). Then one integrates over e 
and e' to obtain the following expression 



drdr' 
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The next step is integration over r and r'. To this end 
we perform the following variables change t = u + v 
and t' = u ~ V with the Jacobian of the transformation 
equaling 2. Denoting y = xe~'^^*e~^'^~^'^ integration 
over u is exact and gives the i5 function S{v — Sihfit). 
Therefore integration over v is also exact. Finally, the 
expectation value reads 



(x^it)) = 2e^ 



(16) 

Here we also use the symmetrical property of the wave 
function. Substituting Eq. ([2]) in Eq. ((T6| . one ob- 
tains that at times t = ti this expression diverges. These 
are the same singularities obtained above in Eq. ([7]). 
Moreover, any "good" Gaussian and exponential func- 
tions lead to these kind of singular behavior for the ex- 
pectation values of physical operators. 

Obviously, these singularities result from a specific 
preparation of the initial wave packets ^'o(a;). Let us pre- 
pare the initial conditions "properly" to obtain the finite 
moments of the physical variables. Owing to Eq. (|12p . 
we present the initial condition as the spectral decompo- 
sition with the Gaussian weight q{e) = [^] * exp(— ae^), 
where real a > 0. This yields the initial wave packet in 
the form of a log-normal distribution 



1 
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with N = Aay/ir. Using Eqs. (fT7|) and (fT6|) . one obtains 
for the nth moment of x 



(x"(i)) 



2ay/TT 
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X / exp [ — ^— J In^ a; + ^ilnx ] d(lna;) 
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exp ( 4wi + a'^n^ + —h^^^t 



(18) 



This behavior of the expectation values is finite and 
spreads exponentially for the arbitrary long time scale. 
This exponential increasing with time consists of two 
values. The first one is the above mentioned classical 



term e 



which is due to the classical motion near 



term is due to the quantum dynamics near the hyper- 
bolic point. 

Both singular behavior and this, completely new, expo- 
nential quantum growrth are due to the nonlinear quan- 
tum dispersion term k = h^T , where T is a characteris- 
tic time scale. For example, for the explosion behavior 
of {SP'it)) in Eqs. and © it is T = For the 

quantum expansion we can define this time scale param- 
eter as h^/ui to define the dimensionless growth of the 



the hyperbolic point. The second, pure quantum, term expectation values exp 



||fi^/i^<^ is dominant and relates to the action of the evo- 
lution operator, which is a delation operator e^^^ , where 
e'"''Pf{x) = f (e-'^^^ (see e.g., 0). Therefore, this 



Lot 



This research was supported by the Israel Science 
Foundation. 



[1] R.S. MacKay in Quantum Chaos, Course CXIX, edited 
by G. Casati, I. Guarneri, U. Smilansky (North-Holland, 
Amsterdam, 1993). 

[2] G.P. Barman, G.M. Zaslavsky, Physica A 91, 450 (1978). 

[3] M.V. Berry, N. Balasz, J. Phys. A: Math. Gen. 12, 625 
(1979) 

[4] M.C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics (Springer- Verlag, Berlin, 1990). 

[5] G. Barton, Ann. Phys. (N.Y.) 166, 322 (1986). 

[6] J. Wang, I. Guarneri, G. Casati, and J. Gong, Phys. Rev. 
Lett. 107, 234104 (2011). 

[7] M.V. Berry, J. P. Keating, in Supersymmetry and Trace 
Formulae, eds. I.V. Lerner, J. P. Keating, D.E. Khmel- 
nitskii, NATO AST Series B 370 (Kluwer Academic, New 
York, 1999). 

[8] J.V. Armitage, in Number Theory and Dynamical Sys- 
tems, eds. M.M. Dodson, J.A.G.Vickers (University Press 
Cambridge 1989). 

[9] R.K. Bhaduri, A. Khare, J. Law, Phys. Rev. E 52, 486 
(1995). 

[10] S. Nonnemacher, A. Voros, J. Phys. A; Mat. Gen. 30, 
295 (1997). 

[11] G.P. Berman, M. Vishik, Phys. Lett. A 319, 352 (2003). 

[12] R.J. Glauber, Phys. Rev. 131, 2766 (1963). 

[13] E. Janke, F. Emde, F. Losh, Tables of Higher Functions 
(McGraw-Hill, NY, 1960). 

[14] The interaction range xo is always a finite value due to 
the local consideration near the hyperbolic point. For ex- 
ample, it can be evaluated as the diiTerence between nor- 
malization constants of the wave functions for the free 
motion Nf and for the saddle point motion Ns (see. Prob- 
lem 133.1 in n^). For the case of free motion the nor- 
malization constant Nf is determined from the integral 

[ii,[a 



/Zo e-'''"/''e''''"/''di = 2'kVTiS{p - p') = Nf5{p - p') . 
Therefore, N; — . We also used here that x has 

the dimension fi^''^. In order to find the A'^s, one uses the 
eigenfunction of the Hamiltonian Hs with an extended 
action on the whole x axis from minus to plus infin- 
ity. Eventually we obtain that the interaction range is 
XQ = Ns-Nf =2Txh^/^. 
[15] L.D. Landau, E.M. Lifshitz, Quantum Mechanics (Perg- 
amon, 1977). 

[16] G. Teschl, Mathematical methods in quantum mechan- 
ics: with applications to Schrodinger operators (Ameri- 
can Mathematical Society, Providence, 2009). 

[17] A mathematically rigorous calculation of the normal- 
ization constant for the wave function X^i^) can be 
pressented by follow the presentation in the monograph 
by V.A. Fock Since the operator xp has continuous 
spectrum e, the eigenfunctions xi^^^) = X^i^) s-rs not 
square integrable. Therefore, the normalization condition 
exists not for the eigenfunction but for the "eigendiffer- 
ential" [l^ Ax(e, a;), which reads 

Ax(^,^)=/r^^x(^',^)*'. 

Substituting here Eq. (|10[) . one obtains 

Ax(e, x) = ^=4— exp[i(e + Ae) In \x\] sin . 

This solution is already square integrable and has the 
following normalization form 
limAE^o ^ /^da;|Ax(e,a;)|^ = Vft. 
Carrying out the following variable change z — 
(Ae/2) In ja::j and taking into account that " dz = 

TT , we obtain A'^ = 47r that coincides exactly with the di- 
mensionless normalization constant Nsh^^^ in Eq. (|10|1 . 
[18] V.A. Fock, Foundations of Quantum Mechanics, (in Rus- 
sian) (Nauka, Moscow, 1976). 



